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Abstract 

We study the nonconserved phase ordering dynamics of the d = 2, 3 
random field Ising model, quenched to below the critical temperature. 
Motivated by the puzzling results of previous work in two and three di- 
mensions, reporting a crossover from power-law to logarithmic growth, 
together with super universal behavior of the correlation function, we 
have undertaken a careful investigation of both the domain growth 
law and the autocorrelation function. Our main results are as follows: 
We confirm the crossover to asymptotic logarithmic behavior in the 
growth law, but, at variance with previous findings, the exponent in 
the preasymptotic power law is disorder- dependent, rather than being 
the one of the pure system. Furthermore, we find that the autocorre- 
lation function does not display superuniversal behavior. This restores 
consistency with previous results for the d = 1 system, and fits nicely 
into the unifying scaling scheme we have recently proposed in the 
study of the random bond Ising model. 
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1 Introduction 



Much recent interest in statistical physics has focused on understanding out- 
of-equilibrium phenomena. In this context, of paramount importance are 
slow relaxation phenomena, which primarily occur in glassy systems. An im- 
portant hallmark of slow relaxation is the lack of time-translation-invariance, 
manifested through aging behavior. A similar phenomenology is also ob- 
served in systems without disorder, e.g., ferromagnets quenched below the 
critical point. The behavior of these systems is well understood in terms of 
the domain growth mechanism of slow relaxation [1, 2J. 

The key feature of domain growth, or coarsening, is the unbounded 
growth of the domain size, which entails scaling due to the existence of a 
dominant length scale, and aging as a manifestation of scaling in multiple- 
time observables. The simplicity of this structure is very attractive and is 
expected to be valid beyond the realm of disorder-free phase-separating sys- 
tems, establishing domain growth as a paradigm of slow relaxation. However, 
as is well known, the applicability of domain growth concepts to hard prob- 
lems (such as spin glasses or structural glasses) still remains a debated issue 
[3]. Therefore, it is of considerable interest to study the role of disorder in 
systems where its presence does not compete with phase ordering [3]. 

A class of systems of this type are disordered ferromagnets, where disorder 
coexists with the low-temperature ferromagnetic order. There are different 
ways to introduce disorder in a ferromagnet without inducing frustration. 
This can be achieved through bond or site dilution, by randomizing the ex- 
change interaction strength while keeping it ferromagnetic, or by introducing 
a random external field. These disordered systems have been an active area 
of research for quite some time now. The unifying theme of investigation 
has been disorder-induced changes in the properties of the underlying pure 
systems, with primary interest in the growth law, in the equal-time corre- 
lation function and, more recently, in the two-time autocorrelation function 
and the related response function El [71 El E] . However, despite the many 
experimental and theoretical studies [4j, a number of issues are still open. 
Among these, of primary importance are (i) the nature of the asymptotic 
growth law (power-law vs. logarithmic) and (ii) the existence of superuniver- 
sal behavior of the correlation and response functions. This is the idea that 
scaling functions are robust with respect to disorder, which is expected not 
to change the low temperature properties of the system [10J. The lack of a 
general framework to understand this complex phenomenology has proven a 
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major obstacle to development. Recently, in the context of the random bond 
Ising model (RBIM) [9J, we have shown that the renormalization group (RG) 
picture of crossover phenomena may well serve the purpose. 

In this paper, we extend the RG conceptual framework to the ordering 
dynamics of the random field Ising model (RFIM) [11] . In this system, the 
deep asymptotic regime turns out to be numerically accessible, allowing us to 
make precise statements regarding the growth law and the super universality 
(SU) issue (see subsection 12. ip . Our principal findings are (a) the existence 
of a crossover from power law domain growth (with a disorder- dependent 
exponent) to logarithmic growth, and (b) the absence of SU. Both results 
fit nicely into an RG picture where disorder acts as a relevant perturbation 
with respect to the pure fixed point. This confirms the robustness and the 
general applicability of the approach proposed in Ref. [9]. 

This paper is organized as follows. In Sec. [21 we provide an overview of 
domain growth laws and of our scaling framework for phase ordering dynam- 
ics in disordered systems. In Sec. El we present detailed numerical results for 
ordering in the d = 2 RFttfl These results are interpreted using the scaling 
framework of Sec. [21 Sec. H] is devoted to the presentation of numerical results 
in the d = 3 case. Finally, in Sec. [51 we conclude this paper with a summary 
and discussion. 

2 Domain Growth Laws in Disordered Sys- 
tems 

Dynamical scaling is the most important characteristic of phase ordering 
systems [H [2]. Let us summarise the concept in the simplest case of pure 
systems. As time increases, the typical domain size L(t) grows and becomes 
the dominant length in the problem. Then, all other lengths can be rescaled 
with respect to L(t). For instance, the two-time order-parameter correlation 
function C(r,t,t w ), with t w < t, scales as (121 [13] 

C(r,t,t w ) = G(r/L,L/L w ) (1) 

where L and L w stand for L(t) and L(t w ), respectively. This contains, as a 
special case, the usual scaling of the equal-time (t = t w ) correlation function 

1 As it is explained in subsection 13. H the spin updating rule we use is equivalent to a 
quench to T = 0. Hence, the d = 2 RFIM phase orders even if d = 2 is the lower critical 
dimensionality. 
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C(r,t) = G\{r/L). Further, for r = 0, we have the aging form of the auto- 
correlation function C(t,t w ) = G 2 (L/L W ). The validity of scaling is, by now, 
a well-established fact. A complete picture of an ordering problem requires 
the understanding of the growth law (i.e. how L(t) depends on t) and of the 
scaling function G(x,y). 

A systematic study of the growth law has been undertaken by Lai et al. 
(LMV) [HI [15] . who identified four different universality classes of growth 
kinetics. LMV considered the role of several factors in the ordering dy- 
namics, e.g., temperature, conservation laws, dimensionality, order param- 
eter symmetry, lattice structure and disorder. An important distinction is 
made between systems that do not freeze (i.e., without free-energy barriers) 
and those that do freeze (i.e., with barriers) when the quench is made to 
T = 0. To the first category belong pure systems with non-conserved dy- 
namics, whose growth follows the power law L(t) = Dt 1 ^, where z = 2. 
LMV designate these as Class 1 systems. To the second category belong 
systems whose growth requires thermal activation. This includes pure sys- 
tems with conserved order parameter and systems (both conserved and non- 
conserved) with quenched disorder. This category is further subdivided into 
three classes. In Class 2 systems the freezing involves only local defects, with 
activation energy Eb independent of the domain size. In this case, growth 
is still power law: L(t) = Dt 1 ^ with z = 2 and 3 for the non-conserved and 
conserved cases, respectively. Furthermore, the prefactor D has a strong tem- 
perature dependence, D ~ e~ EB ^ zT \ Finally, in Class 3 and Class 4 systems, 
the freezing involves a collective behavior which depends on the domain size 
L. If the corresponding activation energy scales with L like Eb{L) ~ eL 93 , 
where e measures the disorder strength, the asymptotic growth law is loga- 
rithmic 

~ (T/e) 1/tp [ln(t/r)] 1/v5 (2) 

with t ~ T/(ipe). For Class 3 systems, we have <p = 1, and for Class 4 
systems, we have (p ^ 1. 

Ferromagnets (with or without disorder) offer examples of the classes 
listed above. For simplicity, let us consider systems with non-conserved order 
parameter. The pure ferromagnetic Ising model with Glauber kinetics is a 
well-known Class 1 system [1J. The d = 1 ferromagnetic RBIM [7] is an 
example of a Class 2 system. The d = 1 RFIM [16] belongs to Class 4, with 
<p = l/2. The RFIM in higher dimensions, d = 2[Hl[T8]andrf = 3[T9l[20], 
shows logarithmic growth, although it is not easy to unambiguously establish 
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the value of (p. Recently, we have presented evidence [5] for logarithmic 
growth in the d = 2 RBIM, but have not established whether it is a Class 3 or 
Class 4 system. This is a particularly interesting system, because its growth 
law was previously |2T] believed to be power law with a disorder-dependent 
exponent. If so, this would have shown the existence of a new universality 
class, say Class 5, in addition to the four listed by LMV. We should stress 
that a huge numerical effort is involved in accessing the logarithmic growth 
regime of the d = 2 RBIM, and our understanding of this system remains 
incomplete. 

In Ref. j9], we have proposed to unify this wide variety of behaviors for 
disordered domain growth into a scaling framework for the growth law itself. 
In all the cases we consider in this paper, disorder (ho) and temperature 
(T) enter through their ratio h /T (see Sec. 13.11 below). This will be de- 
noted by e and, for short, will be termed as disorder. Let us begin with the 
straightforward crossover set-up, where the growth law is assumed to scale 
as 

L{t,e)=&'F(e/1*), (3) 

z = 2 is the growth exponent for non-conserved dynamics in a pure ferro- 
magnet and (f) is the crossover exponent. With the additional assumption 
that the scaling function behaves as 

, , J const., for x — > 0, , , 

iVW^-V*), for x -> oo, { ' 

where x — e/t^, Eq. (j3J) describes the crossover from the power law Lit) ~ t x l z 
to the asymptotic form Lit) ~ t (t/e 1 ^), if <p < 0, and vice-versa if (f) > 0. 
Alternatively, disorder is asymptotically relevant when cf) < 0, and irrelevant 
when <p > 0. The key quantity in the analysis of crossover is the effective 
growth exponent 

1 d\nL(t,e) 1 ,91nF(x) 

(5) 



z e fi (t,e) d\nt z d In x ' 

which depends on t and e through x. 

In the following discussion, it will be useful to use the above relations in 
the inverted form: 

t = L z g(L/X), (6) 

where 

A = e 1/( ^ } (7) 
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is a length scale associated with disorder. The scaling functions appearing 
in Eqs. (j3J) and are related by 

g(y) = F' z {x) (8) 

and y = Lj A is related to x by 

y = x - l '^F{x). (9) 

Then, from Eq. (J3J) and < 0, it follows that 

, v / const., for y < 1, , 
9{v) ~\V-t-\v\ for y»l, (10) 

where i~ x stands for the inverse function of I. The opposite behavior holds 
for (f) > 

\ const., for y > 1. v ; 

Finally, the effective exponent as a function of y is obtained from Eq. fl6]) 

/ \ . d\ng(y) 

ZesiV) = Z + — . 12 

amy 

Therefore, for < (disorder relevant), Eq. fllOp yields 

, s J z, for y < 1, , , 

^efflZ/J | d\D.t-\y)ld\ny, for y > 1, 1 j 

and for > (disorder irrelevant), we obtain from Eq. (ITT]) 

[ for y 3> 1. 
2.1 Superuniversality 

One would expect that the above crossover scenario, which is well estab- 
lished for the growth law, would extend also to the other observables. How- 
ever, this expectation is in conflict with the SU statement that all disorder 
dependence in observables other than the growth law can be eliminated by 
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reparametrization of time through L(t,e) [10] • Thus, according to SU, for 
the autocorrelation function one should have 

C{t, t w , e) = G 2 {L{t w , e)/L{t w , e)) (15) 

where G2 is the scaling function of the pure case. The validity of SU is 
controversial, since the d — 1 results [7J [16] clearly demonstrate the absence 
of SU, while from the study of the correlation function for d > 2, there is 
evidence both in favour [Tj2 1201 Ell E2] and against [0] SU validity Recently, 
the validity of SU has been extended to the geometrical properties of domain 
structures [23] . 

In the next sections we present comprehensive numerical results from 
large scale simulations of ordering dynamics in the RFIM in d = 2, 3. We 
will analyze numerical results within the above scaling framework, producing 
evidence against SU validity. 

3 Numerical Results for d = 2 
3.1 Simulation Details 

We consider an RFIM on a two-dimensional square lattice, with the Hamil- 
tonian 

N 

H =-J} j OiOj - hiGi, a = ±1, (16) 

(ij) <=1 

where (ij) denotes a nearest-neighbour pair, and J > is the ferromagnetic 
exchange coupling. The random field hi = ±/t is an uncorrelated quenched 
variable with a bimodal distribution 

P(hi) = l - [5{h - ho) + 5{h + ho)] . (17) 

The system evolves according to the Glauber kinetics, which models non- 
conserved dynamics [2], with spin flip transition rates given by 

w(<n -> -a,) = i (1 - ^tanh ((flf + ht)/T)) (18) 

where is the local Weiss field. All results in this paper correspond to 
the limit T — > (J/T — > oo), while keeping the ratio e = Hq/T finite. In this 
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limit the system undergoes phase ordering in any dimension, down to d = 1 
[To] . The transition rates take the form 



( 1 for Hfcji < 0, 
w(a t -> -<Ji) = I for Ef ^ > 0, (19) 
{ | (1 - signidihi) tanh(e)) for Hf = 

which shows that disorder affects the evolution through the ratio e = h Q /T, 
as anticipated in Sec. [2j Moreover, Eq. f|T9|) allows for an accelerated up- 
dating rule, with a considerable increase in the speed of computation [24], 
by restricting updates to the sites with H^ai < 0, whose number decreases 
in time as 1/L(t). The gain in the speed of computation becomes more 
important the longer the simulation. 

All statistical quantities presented here have been obtained as an average 
over N Tun = 10 independent runs. For each run, the system has different 
initial condition and random field configuration. We have considered the 
values of disorder amplitude e = 0, 0.25, 0.5, 1, 1.5, 2, 2.5 and we have carefully 
checked that no finite size effects are present up to the final simulation time 
when N = 8000 2 spins. In the pure case, since coarsening is more rapid, we 
have taken N = 12000 2 . 

Numerical results for the growth law and the autocorrelation function are 
presented in the following subsections. 

3.2 Growth Law 

We have obtained the characteristic L(t) from the inverse density of defects. 
This is measured by dividing the number of sites with at least one oppositely- 
aligned neighbor by the total number of sites^]. The plot of L(t, e) vs. t, 
in Fig. [1], shows the existence of at least two time-regimes, separated by 
a microscopic time t of order 1. In the early-time regime (for t < t ), 
there is no dependence on disorder and growth is fast. This is the regime 
where the defects seeded by the random initial condition execute rapid motion 
toward the nearby local minima. For e > and t > to, there is a strong 
dependence on disorder producing slower growth and deviation from the 
power law behavior of the pure case (top cicles line in Fig. [T|). 

2 We have checked that the same results are obtained measuring L(t) from the equal 
time corrlation function. 
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Figure 2: (Color online) Effective exponent z c g vs. t in d = 2. The horizontal 
dashed lines indicate z, the plateau values of z e s. 



In Fig. El we show the time-dependence of the effective exponent z e g(t, e), 
defined by Eq. (jHJ). For t > to, this plot shows the existence of an intermediate 
power-law regime, characterized by a plateau where z e s is approximately 
constant. This is followed by the late regime where z e g is clearly time- 
dependent. The disorder dependent values of z efi on the plateaus, denoted 
by z, are listed in Table [1] and plotted in Fig. [3j We encountered a similar 
crossover in our study of the d = 2 RBIM [9], i.e., a preasymptotic power 
law regime with a disorder dependent exponent, followed by an asymptotic 
regime where the growth law deviates from a power law. 

The appearance of a disorder dependent exponent z in the intermediate 
regime suggests to upgrade the crossover picture, presented in Sec. El by 
replacing the pure growth exponent z by z in all the scaling formulae. Then, 
from Eq. ([121) it follows that z e s—z ought to depend only on y = L/X. Indeed, 
as Fig. H] shows, it is possible to determine numerically the quantity A such 
that the plots of (z e g — z) vs. L/X, for different disorder values, collapse on a 
single master curve. The e-dependence of A is displayed in Fig. [5] and is well 



10 



IN 




0.4 

e 



Figure 3: (Color online) z (taken from Fig. [2]) vs. e. The red line is the best 
fit z = 2.0 + 1.4e L35 . 



e 


z 





2 


0.2 


2.20 


0.35 


2.31 


0.5 


2.52 


0.65 


2.77 


0.8 


3.05 



Table 1: Plateau exponent z for various disorder strengths. 
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Figure 4: (Color online) Subtracted effective exponent (z e g — z) vs. L/X. 
The dashed line is the best fit z e g — z = 0.0055(L/A) 1 ' 5 . 



fitted by 

A ~ e- 2 . (20) 

Comparing this with Eq. ([7]), the negative exponent implies < and, 
therefore, that disorder acts like a relevant scaling field. This is also confirmed 
by the behavior of z eS (y) in Fig. HJ which is consistent with Eq. (TC?j) but not 
with Eq. (HI]). 

Fitting the data of Fig. H]to the power law z cff — z = by 1 ?, we find b ~ 
0.0055 and (p ~ 1.5. Hence, from the definition of z c q in Eq. ()5]) follows 

dint , ,„ . , 

— = z + W (21) 

which, after integrating with respect to L, yields 

t = K(e)I/g{L/\) (22) 
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Figure 5: (Color online) Plot of A vs. e. 



e 


K 


0.5 


0.39 


1.0 


0.13 


1.5 


0.039 


2.0 


0.015 



Table 2: Prefactor K(e) for various disorder strengths. 

where K(e) is an e-dependent prefactor. Indeed, replotting in Fig. [6]the data 
of Fig. [U as tL~ z /K(e) vs. y, an excellent data collapse on the master curve 

g(y) ~ exp ( 23 ) 

is obtained, with the values of K(e) listed in Table HJ 

The plot of the scaling function g{y) illustrates quite effectively (i) the 
existence of the crossover, and (ii) that our numerical data reach deep into 
the asymptotic regime. The flat part of the curve, where g(y) lies on the 
horizontal dashed line at g(y) = 1, corresponds to the preasymptotic power- 
law regime [cf. Eq. ffTUj) ]. The sharp and fast increase of g(y), for large y, 
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corresponds to the crossover to the asymptotic growth law 



L 



|ln(t/A*) 



(24) 



A 



which corresponds to the Class 4 form of Eq. (J2J) . 

Summarising, our main findings for the growth law, in the d = 2, case 
are as follows: 

1. Disorder is a relevant perturbation with respect to pure- like behavior. 

2. The corresponding growth law shows a clear crossover from power- law 
to logarithmic behavior: 



This differs from previously found results, since the preasymptotic 
power law is not pure-like, due to the e-dependence of the exponent 
z. This feature, also observed in the d = 2 RBIM [9], means that dis- 
order although globally relevant acts like a marginal operator in the 
neighborhood of the pure fixed point [25J. 

Finally, we remark on the considerable numerical advanyage in using the 
effective exponent as a probe for the crossover. In fact, while the switch from 
preasymptotic to asymptotic behaviors in z e g takes place at about L cr ~ A, 
from Eqs. (1221 and ([23]) . it follows that the condition by^/ip = 1 puts the 
crossover, when looking at the domain size, at the much greater value L cr ~ 
50 A, as it is evident from Fig. [6j 

3.3 Autocorrelation Function and SU Violation 

The results presented above show that disorder affects the growth law as 
an asymptotically relevant parameter. Therefore, one would expect this to 
apply also to other observables. However, as explained in Sec. 12.11 such an 
expectation would be in conflict with claims of SU validity. 

In this section, we study the autocorrelation function, defined by 




(25) 




(26) 
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which is independent of i, due to space-translation invariance. In Fig. [7] 
C(t,t w ,e) has been plotted against L/L w , for e = 0,0.5,0.65 and with dif- 
ferent values of t w , chosen in such a way that the ratio v = L w /X takes the 
three different values v — 0, 0.25, 0.85. If SU were valid all the curves, ir- 
respective of the value of e, should collapse on the e = 0, or v = 0, master 
curve. Instead, there is an evident e-dependence which excludes SU validity. 
In addition, curves with the same value of v do collapse, showing that the 
autocorrelation function obeys the extended aging form 



The pure master curve, with v = 0, lies below the two corresponding to 
v = 0.25 and v = 0.85. Figure |8] displays a similar plot, with v = and 
v = 0.25. The latter value is obtained by combining the four different values 
of disorder e = 0.2, 0.35, 0.5, 0.65 with appropriately chosen t^-values. Again, 
there are two distinct master curves for different v-values. 

The reported SU violation is in agreement with the behavior of the auto- 
correlation function in the d — 1 RFIM [16] . 

4 Numerical Results for d = 3 

As mentioned above, in previous studies of the d = 3 RFIM [TjJJ 120] SU has 
been found to hold. Here, instead, we present new results for this system, 
which produce evidence for the same pattern of SU violation observed in the 
d = 2 case. 

Simulations were made on a system with N = 300 3 spins on a cubic 
lattice, evolving with the transition rates f[T§j) and averaging over iV run = 20 
runs, each with different initial condition and random field configuration. We 
have considered the three disorder values e = 0.5, 1, 2. Although the quality 
of the data does not allow for an analysis of high precision as in the d = 2 
case discussed above, nonetheless the main features, including the lack of SU 
validity, do emerge quite clearly. 

Let us begin with the growth law. The L(t) data have been plotted in 
Fig. |9j The qualitative behaviour is the same as in Fig. [Tj namely, as 
disorder increases, growth slows down. The corresponding effective exponent 
z e fi is displayed in Fig. [TTJJ For the smaller values, e = 0.5 and e = 1, 
the overall behaviour is qualitatively similar to that of Fig. 121 showing a 




(27) 
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Figure 7: (Color online) Autocorrelation function in d — 2, for disorder 
values and waiting times t w chosen so that v = L w /X takes the three values 
v = 0, 0.25, 0.85 corresponding, from bottom to top, to the three different 
master curves. 
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Figure 9: (Color online) Growth law in d = 3 for different disorder values. 



crossover from power law (with disorder dependent exponent) to logarithmic 
behavior, as in Eq. f )25|) . The data for the highest disorder value e = 2, 
instead, display a novel behaviour. The preasymptotic power law regime 
disappears and is replaced by a pronounced peak in the effective exponent, 
which reveals strong pinning of the interfaces. Therefore, the nature of the 
crossover is qualitatively different for weak and for strong disorder. The 
detailed investigation of this novel feature is delayed to a future publication. 

The data for the autocorrelation function are displayed in Fig. [TTJ As 
stated above, the quality of the data is not sufficient to carry out a precise 
scaling analysis as in Sec. El In particular, it is not possible to extract 
the characteristic length A reliably and to organise the plot with a choice 
of t w values aimed to keep constant the ratio L w /X. However, the evident 
e-dependence in Fig. [11] is quite sufficient to exclude SU validity. 

As a final remark, it should be noted that with the algorithm illustrated 
in Sec. 13.11 we consider a quench to T = 0, while in Refs. [HJ [20] quenches 
to T > were considered. Therefore, it remains an open question, to be 
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Figure 10: (Color online) Effective exponent z e s{t,e) in d = 3. 
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Figure 11: (Color online) Autocorrelation function in d = 3 for different e 
and t w values. 



investigated, whether the discrepancy between ours and previous results may 
be related to the role of the final temperature. 

5 Summary and Discussion 

Let us conclude this paper with a summary and discussion of our results. We 
have recently initiated a large-scale simulation study of nonconserved domain 
growth in disordered systems. Our study was motivated by some ambiguities 
existing in the available literature: 

(a) The precise nature of the asymptotic growth law in the standard models 
(RBIM, RFIM, etc.) was unclear. 

(b) Quantities like the equal-time correlation function (or its Fourier trans- 
form, the structure factor) showed SU, i.e., the scaling functions were inde- 
pendent of disorder. It was not clear to us why the crossover in the domain 



21 



growth law was not accompanied by a corresponding crossover in the corre- 
lation function. 

(c) There were very few studies of two-time quantities like the autocorrelation 
function and the response function. 

With this background, we investigated two-time quantities in the d — 1, 2 
RBIM [3 E] and, with the present paper, we have extended the investigation 
to the d — 2, 3 RFIM with Glauber spin-flip kinetics. Our results can be 
summarized as follows: 

(i) First, we have formulated a general scaling framework for the study of 
disordered domain growth. The framework is based on the RG concept of 
asymptotically relevant parameter, and proved very convenient for interpret- 
ing our earlier RBIM results [H] . In this paper, we have used this framework 
to successfully understand crossovers in the RFIM. 

(ii) Second, we find that there is a crossover in the domain growth law from a 
preasymptotic regime showing power law growth with a disorder dependent 
exponent \L{t) ~ t 1 ^^] to an asymptotic regime with logarithmic growth 
[L(t) ~ (lnt) 1 /* 3 with (p ~ 1.5]. Following the analysis of Ref. [2T], this can 
be related to an underlying crossover from logarithmic dependence of the free 
energy barriers on the domain size to power law dependence. The mechanism 
producing the power law dependence is particularly clear in d = 1, where the 
interface motion can be mapped into the random walk in a random potential 
of the Sinai type pHETj. 

(iii) Third, and perhaps most important, we find that the autocorrelation 
function does not obey SU. The scaling function shows a crossover corre- 
sponding to the crossover in the growth law. 

In the light of the above results, what is the path ahead? We are currently 
investigating other disordered systems to confirm whether the domain growth 
scenario is consistent with the scaling picture developed here. It is also rele- 
vant to re-examine earlier results demonstrating SU for the equal-time corre- 
lation function. It is possible that the equal-time correlation function has a 
delayed crossover, and may violate SU in larger and longer simulations. Al- 
ternatively, it could be that the equal-time correlation function is a relatively 
crude and feature-less characteristic of the morphology. It may be worthwhile 
to study more sophisticated measures of the morphology [26], which could 
show differences between pure and disordered phase ordering systems. Our 
studies demonstrate that there remain many unanswered questions in this 
area. We hope that our work will motivate fresh interest in these problems. 
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